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Abstract. We first prove that if a has a prime factor not dividing b then there are 
infinitely many positive integers n such that ( Q "+ b ™) is not divisible by bra + 1. This 
confirms a recent conjecture of Z.-W. Sun. Moreover, we provide some new divisibility 
properties of binomial coefficients: for example, we prove that (3^) and (^*) are divis- 
> ible by 6ra — 1, and that (^g ^ 1 ) is divisible by 66ra — 1, for all positive integers n. As we 

\& ! show, the latter results are in fact consequences of divisibility and positivity results for 

quotients of q-binomial coefficients by ^-integers, generalizing the positivity of g-Catalan 
numbers. We also put forward several related conjectures. 
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I. Introduction 

The study of arithmetic properties of binomial coefficients has a long history. In 1819, 
d • Babbage [6] proved the congruence 

2p- T 
p — 1 

for primes p ^ 3. In 1862, Wolstenholme [28] showed that the above congruence holds 
modulo p 3 for any prime p ^ 5. See [20] for a historical survey on Wolstenholme's theorem. 
Another famous congruence is 

'2n\ 

= (mod n + 1). 



n 
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The corresponding quotients, the numbers C n := ^-j- ( 2 ™) , are called Catalan numbers, and 
they have many interesting combinatorial interpretations; see, for example, [12] and [24, 
pp. 219-229]. Recently, Ulas and Schinzel [27] studied divisibility problems of Erdos and 
Straus, and of Erdos and Graham. In [25, 26], Sun gave some new divisibility properties 
of binomial coefficients and their products. For example, Sun proved the following result. 

Theorem 1.1. [26, Theorem 1.1] Let a, b, andn be positive integers. Then 

an + bn\ bn + 1 

= mod — - — ; — . (1.1) 

an J gcd(a, bn+1) 

Sun also proposed the following conjecture. 

Conjecture 1.2. [26, Conjecture 1.1] Let a and b be positive integers. If 
(bn + 1) \ ( an ^ >n ) for all sufficiently large positive integers n, then each prime factor of 
a divides b. In other words, if a has a prime factor not dividing b, then there are infinitely 
many positive integers n such that (bn + 1) \ { an ^ n ) • 

Inspired by Conjecture 1.2, Sun [26] introduced a new function /: Z + x Z + — > N. 
Namely, for positive integers a and b, if ( an a + bn ) is divisible by bn + 1 for all n e Z+, then 
he defined f(a, b) = 0; otherwise, he let f(a, b) be the smallest positive integer n such that 
( a? on ™) * s n °t divisible by bn + 1 . Using Mathematica, Sun [26] computed some values of 
the function /: 

/(7, 36) = 279, /(10, 192) = 362, /(ll, 100) = 1187, /(22, 200) = 6462, .... 

The present paper serves several purposes: first of all, we give a proof of Conjecture 1.2 
(see Theorem 2.1 below); second, we provide congruences and divisibility results similar 
to the ones addressed in Theorem 1.1 and Conjecture 1.2 (see Theorems 2.2-2.3 in Sec- 
tion 2); third, we show in Section 3 that among these results there is a significant number 
which can be "lifted to the q- world;" in other words, there are several such results which 
follow directly from stronger divisibility results for g-polynomials. In particular, Theo- 
rem 1.1 is an easy consequence of Theorem 3.3, and Theorem 2.3 is an easy consequence 
of Theorem 3.1. On the other hand, Theorem 2.4 hints at the limitations of occurrence 
of these divisibility phenomena. Sections 4-6 are devoted to the proofs of our results in 
Sections 2 and 3. We close our paper with Section 7 by posing several open problems. 

2. Results, I 

Our first result is a more precise version of Conjecture 1.2. 

Theorem 2.1. Conjecture 1.2 is true. Moreover, if p is a prime such that p \ a butp\b, 
then 

a + b 

where <p(n) is Euler's totient function. 

For the proof of the above result, we need the following theorem. 
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Theorem 2.2. Let a and b be positive integers with a > b, and f3 an integer. Let p be a 
prime not dividing a. Then there are infinitely many positive integers n such that 



an 
bn + (} 



±1 (modp). 



Our proofs of Theorems 2.1 and 2.2 are based on Euler's totient theorem and Lu- 
cas' classical theorem on the congruence behaviour of binomial coefficients modulo prime 
numbers, see Section 4. 



In [14, Corollary 2.3], the first author proved that 

eO (mod2n-l). (2.1) 



6n 
3n 



It is easy to see that 



(2n\ (2n - 1\ 4n - 2 /2n - 2\ . . . . 

= 2 = ,1=0 mod2n-l). 2.2 

\n J \ n J n \ n — 1 J 



The next theorem gives congruences similar to (2.1) and (2.2). 
Theorem 2.3. Let n be a positive integer. Then 



12n\ _ (Yin 
3n J \4n 

30n 
5n 



= (mod 6n-l), (2.3) 



eO (mod (lOn- l)(15n- 1)), (2.4) 
Q0n\ _ (120n\ _ /120n\ 



. _ . , - , . = (mod 30n- 1), (2.5) 



/330n\ 
\88n ) 



(mod 66n-l). (2.6) 



We shall see that this theorem is the consequence of a stronger result for g-binomial 
coefficients, cf. Theorem 3.1 in the next section. 

It seems that there should exist many more congruences like (2.1)-(2.6). On the other 
hand, we have the following negative result. 

Theorem 2.4. There are no positive integers a and b such that 

( m + bn )=0 (mod3n-l) 
\ an J 

for all n ^ 1 . 

For a possible generalisation of this theorem see Conjecture 7.2 in the last section. 
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3. Results, II: ^-divisibility properties 
Recall that the q-binomial coefficients (also called Gaufiian polynomials) are defined by 

(l-g n )(l-g n_1 )---(l-g) 



q k )(l-q 



fc-D 



;i_g)(l_gn-fc)(l_g. 



if ^ k ^ n, 
otherwise. 



We begin with the announced strengthening of Theorem 2.3. 
Theorem 3.1. Let n be a positive integer. Then all of 



1 



Q 



7 6n-l 



12n 
3n 



1 - a 6 ™" 1 

i-q 

I _ g30n-l 



Yin 
An 

"I20n 
40n 



1-q 



7 30n-l 



1 - q° 
1-g 

' X _ g30n-l 



60n 
6n 



i 

120n 
45n 



1-q 



1-q 



66n-l 



330n 
88n 



are polynomials in q with non-negative integer coefficients. Furthermore, 



(1 -glOn-l)(l_g 



15n-l^ 



30n 
5n 



(3.1) 



(3.2) 



is a polynomial in q. 



It is obvious that, when a = b = 1, the numbers ( an ^ >n ) /(bn + 1) (featured implicitly in 
Conjecture 1.2 and in Theorem 2.1) reduce to the Catalan numbers C n . There are various 
g-analogues of the Catalan numbers. See Furlinger and Hofbauer [10] for a survey, and 
see [11, 17, 16] for the so-called q, t-Catalan numbers. 

A natural g-analogue of C n is 



C n (q) 



1 



1 - q n + x 



In 
n 



It is well known that the g-Catalan numbers C n (q) are polynomials with non-negative 
integer coefficients (see [1, 2, 4, 10]). Furthermore, Haiman [17, (1.7)] proved (and it 
follows from Lemma 5.2 below) that the polynomial 



-,bn+l 



bn + n 
n 



has non-negative coefficients for all b, n ^ 1. Another generalisation of C n (q) was intro- 
duced by the first author and Zeng [15]: 



1 - q n 



In 




r 2n - f 




In -I 


n — k 




n — k 


i 


n — k — 1 



q\ 



1 < k < n. 



They noted that the -B n; fc(g)'s are polynomials in q, but did not address the question 
whether they are polynomials with non-negative coefficients. As the next theorem shows, 
this turns out to be the case. The theorem establishes in fact a stronger non-negativity 
property. 
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Theorem 3.2. Let n and k be non-negative integers with ^ k ^ n. Then 



1-q 



gcd(fc,n) 



Q 



2n 
n — k 



(3.3) 



is a polynomial in q with non-negative integer coefficients. In particular, B nt k(q) = 



l-g k r I 
l-q n in 



2n 



k \ q is a polynomial with non-negative coefficients. 
Applying the inequality [26, (2.1)], we can also easily deduce that 



a 



a,b,n{q) '■— ^ _ 



l-q° 



bn+1 



an + bn 
an 



is a product of certain cyclotomic polynomials, and therefore a polynomial in q. Again, 
as it turns out, all coefficients in these polynomials are non- negative. Also here, we 
have actually a stronger result, given in the theorem below. It should be noted that it 
generalises Theorem 1.1, the latter being obtained upon letting q — > 1. 



Theorem 3.3. Let a, b, and n be positive integers. Then 



q 



gcd(an,6n+l) 



-,bn+l 



an + bn 
an 



q 



gcd(an,6n+l) 



1 an+bn+l 



an + bn + 1 
an 



is a polynomial in q with non-negative coefficients. 
Corollary 3.4. Let a, b, and n be positive integers. Then 



l-q a 



an + bn 
an 



I _ gbn+l 

is a polynomial in q with non-negative coefficients. 

The proofs of the results in this section are given in Section 5. 



4. Proofs of Theorems 2.1 and 2.2 

The proof of Theorem 2.2 (from which subsequently Theorem 2.1 is derived) makes 
essential use of Lucas' classical theorem on binomial coefficient congruences (see, for 
example, [7, 9, 13, 21]). For the convenience of the reader, we recall the theorem below. 

Theorem 4.1 (LUCAS' THEOREM). Let p be a prime, and let a , b , ■ ■ ■ , a m , b m G {0, 1, 
...,p-l}. Then 

fa + a x p H h a rn p m \ -pr f Oi\ , , N 

U+m.+-+^J s SU>) ( p) - 

Proof of Theorem 2.2. Note that gcd(p, a) = 1. By Euler's totient theorem (see [23]), we 
have 

p<p(o) _ ! = o (mod a). 
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Since a > b > 0, there exists a positive integer N such that an > bn + (3 > holds for all 
n > N. Let r be a positive integer such that p ri P( a ) — 1 > aN, and let n = (jf^i -) — I) /a. 
Then, by Lucas' theorem, we have 

where bn + /3 = b + b\p + ■ ■ • b m p m with ^ b 0: . . . , b m ^ p — 1. It is clear that there are 
infinitely many such r and n. This completes the proof. □ 

Proof of Theorem 2.1. Suppose that a and b are positive integers and p a prime such that 
p | a but p\b. We have the decomposition 

an + on\ f an + bn\ a + b fan + bn — 1\ . , . 

(4.1) 



6n + l\ an J \ an — 1 / a \ an — 2 

It is clear thatpf (a+6). By the proof of Theorem 2.2, if we take n = (p rip(a+b) -l)/(a + b) 
(r ^ 1), then 

'an + bn\ 

±1 (mod p), 



an — 1 



and thus 



Z 



, s /an + bn — 1\ an — 1 fan + 6rA . , . , . , . 

a + 6 = =±a + 6#0 mod p. 4.2 

\ an - 2 / n \ an — 1 J 

Combining (4.1) and (4.2) gives 

1 fan + bn" 
bn + 1 \ an 

for all n = (p r ^ a+b "> - I) /{a + 6) (r = 1, 2, . . .). Namely, Conjecture 1.2 holds and 

/(a,6)^^-— — , 
a + o 

as desired. □ 



5. Proofs of Theorems 3.1-3.3 and of Corollary 3.4 

All the proofs in this section are similar in spirit. They all draw on a lemma from [22, 
Proposition lO.l.(iii)], which extracts the essentials out of an argument of Andrews [3, 
Proof of Theorem 2]. (To be precise, Lemma 5.1 below is a slight generalisation of [22, 
Proposition lO.l.(iii)]. However, the proof from [22] works also for this generalisation. We 
provide it here for the sake of completeness.) Recall that a polynomial P(q) = Yl^oPiQ 1 
in q of degree d is called reciprocal if Pi = Pd-i for all i, and that it is called unimodal if 
there is an integer r with ^ r ^ d and ^ po ^ . . . ^ p r ^ . . . ^ pd ^ 0. 

Lemma 5.1. Let P(q) be a reciprocal and unimodal polynomial and m and n positive 
integers with m ^ n. Furthermore, assume that A(q) = \ q qn P(q) is a polynomial in q. 
Then A(q) has non-negative coefficients. 
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Proof. Since P(q) is unimodal, the coefficient of z k in (f — q m )P(q) is non-negative for 
^ k < deg(P)/2. Consequently, the same must be true for A(q) = j^P(q), considered 
as a formal power series in q. However, also A(q) is reciprocal, and its degree is at most the 
degree of P(q). Therefore the remaining coefficients of A(q) must also be non-negative. □ 

Proof of Theorem 3.1. In view of Lemma 5.1 and the well-known reciprocality and uni- 
modality of g-binomial coefficients (cf. [24, Ex. 7.75 .d]), for proving Theorem 3.1 it suffices 
to show that the expressions in (3.1) and (3.2) are polynomials in q. We are going to ac- 
complish this by a count of the cyclotomic polynomials which divide numerators and 
denominators of these expressions, respectively. 



We begin by showing that j 
fact that 



_ n 6n — 1 



12n] 
3n \q 



is a polynomial in q. We recall the well-known 



l-q n = H^ d (q), 

d\n 

where §d(q) denotes the d-th cyclotomic polynomial in q. Consequently, 



6n-l 



l-q 



Yin 
3n 



Yin 



d=2 



with 



12n 




3n 




9n 


d 




d 




d 



(5.1) 



ed = ~x(d I (6n - 1)) + 

where x(S) = 1 if 5 is true and x($) = otherwise. This is clearly non-negative, unless 
d | (6n-l). 

So, let us assume that d \ (6n — 1), which in particular means that d ^ 5. Let us write 
X = {3n/d}, where {a} := a — [a\ denotes the fractional part of a. Using this notation, 
Equation (5.1) becomes 

e d = - X (d | (IdX - 1)) + L4XJ - L3XJ . 

Since ^ X < 1, we have — 1 ^ IdX — 1 < Id. The only integers which are divisible 
by d in the range —1, 0, . . . , 2d — 1 are and d. Hence, we must have X = 1 / (2c?) or 
X = [d + 1)/ {2d). The former is impossible since X is a rational number which can be 
written with denominator d. Thus, the only possibility left is X = (d + I)/ (2d). For 
this choice, it follows that |_4XJ — |_3XJ =2 — 1 = 1. (Here we used that d > 5.) This 
proves that is non-negative also in this case, and completes the proof of polynomiality 

of l ~1 \ 12n 1 
Ui i_ g 6n-l [ 3n \q- 

The proof of polynomiality of 1 _ 1 g ^n-i [ 1 in] q * s completely analogous and therefore left 
to the reader. 

We next turn our attention to [ 6 6 °^] a - Again, we write 



l-q 



1 — q 



30n-l 



60n 
6n 



60n 



d=2 
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with 



e d = -x{d\ (30n-l)) + 



QOn 




6n 




54n 


d 




d 




d 



(5.2) 



This is clearly non-negative, unless d | (30n — 1). 

We assume d | (30n — 1) and note that this implies d = 7 or d ^ 11. Here, we write 
X = {6n/d}. Using this notation, Equation (5.2) becomes 

e d = ~x(d I (5dX - 1)) + L10XJ - L9XJ . 

Since ^ X < 1, we have d \ (5dX — 1) if and only if X is one of 

11 12 13 14 1 

5c?' 5 5c?' 5 5cf 5 5cf 5 5d 

For the same reason as before, the option X = l/(5d) is impossible. For the other options, 
the corresponding value of [lOJ^J — |_9J^J is always 1, except if X = | + ^ and d = 7. 
However, in that case, we have X — Jjr, which cannot be written with denominator 
d = 7. Therefore this case can actually not occur. This completes the proof that e d is 
non-negative for all d ^ 2, and, hence, that [ 6 6 °™] 9 is a polynomial in q. 

Proceeding in the same style, for the proof of polynomiality of [ 
to show that 

e d = - x {d | (3dX - 1)) + [12XJ - [4XJ - [8XJ (5.3) 

is non-negative for all X = x/d with ^ x < d, x being integral. Clearly, the expression in 
(5.3) is non-negative, except possibly if d \ (3dX — 1). With the same reasoning as before, 
we see that the only cases to be examined are X — (d+ l)/(3c?) and X = (2d + l)/(3c?), 
where d = 7 or d ^ 11. If X = (d + l)/(3d), we have 



40n \ g we have 



[12XJ - [4XJ - [8XJ =4-1 



3 + 3d 



So, this will be equal to 1, except if d — 7. However, in that case we have X — 
which cannot be written with denominator d = 7, a contradiction. Similarly, if X = 
(2d+ l)/(3d), we have 

[12XJ - [4XJ - [8XJ = 8- 2- 5 = 1. 

So, again, the exponent in (5.3) is non-negative, which establishes that 1 _^n-i l^on ] q 
is a polynomial in q. 

For the proof of polynomiality of [ 3 88°n ] q we have to show that 

e d = ~x(d | {3dX - 1)) + [15XJ - [4XJ - [11XJ (5.4) 

is non-negative for all X = x/d with ^ x < d, x being integral. Clearly, the expression in 
(5.4) is non- negative, except possibly if d | (3dX — 1). With the same reasoning as before, 
we see that the only cases to be examined are X = (d + 1)/ (3d) and X = (2d +1)/ (3d), 
where d = 5, d = 7, or d ^ 11. If X = (d+ l)/(3d), we have 

L15XJ - [4XJ - = 5 + 



5 


- 1 - 


n 


ii 

+ 3d_ 


A 




y 
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So, this will be equal to 1, except if d — 7. However, again, this is an impossible case. 
Similarly, if X = (2d + I) /(3d), we have 



So, again, the exponent in (5.4) is non-negative, which establishes that [ 3 8 3 8 °™ ] q 



[15X\ - [AX\ - [11X\ =10 + 

expon< 
is a polynomial in q. 

Turning to (3.2), to prove that (^lon-^^isn-i) 
show that 



22 11 
T + 3d 



= 1. 



3 5n] q 1S a polynomial in q, we must 



e d = - x [d | (2dX - 1)) - x (d | (3dX - 1)) + [6XJ - [5XJ 



(5.5) 



is non-negative for all X = x/d with ^ x < d, x being integral, d ^ 5. 

First of all, we should observe that gcd(10n — 1, 15n — 1) = 1, whence the two truth 
functions in (5.5) cannot equal 1 simultaneously. Therefore, the expression in (5.5) is non- 
negative, except possibly if d \ (2dX — 1) or if d \ (3dX — 1). With the same reasoning as 
before, we see that the only cases to be examined are X = (d + 1)/ (2d), X = (d+1)/ (3d), 
and X = (2d + l)/(3d), where, in the latter two cases, the choice of d = 3 is excluded. If 
X = (d+1)/ (2d), then 



3 




5 

.2 + 


5 


A 




2d_ 



[6X\ - [5X\ =3 + 
which always equals 1 for d ^ 2. If X — (d+ I) /(3d), then 

L6XJ - [5XJ = 2 + 



2 




5 

j + 


5 


A 




3d_ 



which always equals 1 for d — 2, 6, 7, . . . (sic!). Because of our assumptions, we do not 
need to consider the cases d = 3 and d = 5, so the only remaining case is d = 4. However, 
in that case X — which cannot be written with denominator d = 4, a contradiction. 
Finally, if X = (2d + I) /(3d), then 

L6XJ - [5XJ = 4 + 

which always equals 1 for d ^ 2. 

So, again, the exponent in (5.5) is non-negative in all cases, which establishes that 



2 




10 


5 


A 




J" 


" 3d_ 



(i-q) 2 



'30ral 



(l-q 10 "- 1 )(l-q 15 "- 1 ) il5nlq 



is a polynomial in q. 



□ 



Proof of Theorem 3.2. By Lemma 5.1, it suffices to establish polynomiality of (3.3). When 
written in terms of cyclotomic polynomials, Expression (3.3) reads 



1 



gcd(fc,n) 



with 



2n 
n — k 



e d = x{d | gcd(fc, n)) - x(d \ n) + 



2n 



i[M<i) ed , 



d=2 



2n 




n — k 




n + k 


. d . 




d 




d 



(5.6) 
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Similar to before, let us write N = {n/d} and K = {k/d}. Using this notation, Equa- 
tion (5.6) becomes 

e d = X (d | gcd(A;, n)) - X (d \ n) + L2iVj - [N - K\ - [N + K\ . (5.7) 

We have to distinguish several cases. If d \ n, then N = 0, and (5.7) becomes 

e d = X (d\k)-l- [-K\ . 

We see that this is zero (and, hence, non- negative) regardless whether d | k or not. 
On the other hand, if we assume that d\n, then (5.7) becomes 

e d = [2N\ - [N - K\ - [N + K\ , 

and this is always non-negative. We have proven that (3.3) is indeed a polynomial in q. 

The statement on B njk (q) follows immediately from the fact that gcd(k,n) \ k. □ 

Finally, Theorem 3.3 follows immediately from the following strengthening of a non- 
negativity result of Andrews [3, Theorem 2]. 



Lemma 5.2. Let a and b be positive integers. Then 



q 



gcd(a,b) 



1 — q 



a+b 



a + b 
a 



(5.8) 



is a polynomial in q with non-negative integer coefficients. 



Proof. In view of Lemma 5.1, it suffices to show that the expression in (5.8) is a polynomial 
in q. Again, we start with the factorisation 



I _ qgcd(a,b) 

1 - q a+b 



a + b 
a 



a+b-l 



n ^ ed > 



d=2 



with 



ed = x(d | gcd(a,6)) + 



a + b-l 




a 




b 


d 




.d- 




d_ 



(5.9) 



Next we write A = {a/d} and B = {b/d}. Using this notation, Equation (5.9) becomes 

1 



ed = x(d | gcd(a,6)) + 



A + B- 



d 



This is clearly non-negative, unless A = B = 0. However, in that case we have d \ a and 
d | b, that is, d \ gcd(a, b), so that is non-negative also in this case. □ 



Proof of Corollary 3.4- This follows immediately from Theorem 3.3 and the fact that 
a | gcd(a, bn + 1) — gcd(an, bn + 1). □ 
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6. Proof of Theorem 2.4 



The following auxiliary result on the occurrence of prime numbers congruent to 2 mod- 
ulo 3 will be crucial. 

Lemma 6.1. If x ^ 3761, there is always at least one prime number congruent to 2 
modulo 3 contained in the interval (x, f§^)- 

Proof. Let 9(x; 3, 2) denote the classical Chebyshev function, defined by 

0(x;3,2)= l °ZP- 

p prime, p^x 
p=2 (mod 3) 

McCurley proved the following estimates for this function (see [19, Theorems 5.1 and 
5.3]): 

%;3,2)<0.51y, y>0, 

9(y; 3, 2) > 0.49 y, y > 3761. 

This implies that, for x > 3761, we have 

9 (§x; 3, 2) - 9(x; 3, 2) > 0.49 • §x - 0.51 x > 0.0057 a; > 1. 

This means that there must be a prime number congruent to 2 modulo 3 strictly between 
x and j^x, which is what we wanted to prove. (To be completely accurate: the above 
argument only shows that such a prime number exists in the half-open interval (x, 
However, existence in the open interval (x, ^x) can be easily established in the same 
manner, by slightly lowering the value of yjj in the above argument.) □ 

Proof of Theorem 2.4- We first verified the claim for a, b ^ 13000 using a computer. 



To establish the claim for the remaining values of a and b, we have to distinguish several 
cases, depending on the congruence classes of a and b modulo 3 and the relative sizes of 
a and b. 

First let (a, b) e {(0,0), (0,1), (1,0), (0,2), (1,1), (2, 0)} + (3Z) 2 . By Dirichlet's theo- 
rem [8] (see [5]), we know that there are infinitely many primes congruent to 2 modulo 3. 
Let us take such a prime p with p > a + b, and let 3n — 1 = p, that is, n — (p + l)/3. 
Furthermore, let v p (a) denote the p-adic valuation of a, that is, the maximal exponent e 
such that p e divides a. Writing a = 3ai + 02 and b = 3&i + b 2 with ^ a 2 , b 2 ^ 2, by the 
well-known formula of Legendre [18, p. 10] for the p-adic valuation of factorials, we then 
have 



3n- 1 



an + bn 
an 



. a 2 + b 2 



a + b)n 



P 



an 



pf 



bn 

pi 



+ 



a + b 



Since, in the current case, we have $C a 2 



3p 

b 2 ^ 2 and 



a 2 a 
~3 + 3p 



3 



b 

3p 



(6.1) 



^ < I, all the values of the floor 

3p 3 ' 

functions on the right-hand side of the above equation are zero. Consequently, the p-adic 
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valuation of 3^zr( an a ^ n ) equals —1 for our choice of p and n, which in particular means 
that ^— ( an+bn ) i s not an integer. 

in— 1 V an / 

For the next case, consider a pair (a, b) e {(2, 1), (2,2)} + (3Z) 2 . Let us first assume 
that b ^ ^a. Since we have already verified the claim for a,b ^ 13000, we may assume 
a ^ 13000. We now choose a prime p = 2 strictly between and a. Such a prime 



is guaranteed to exist by Lemma 6.1, because, due to our assumption, we have ^a ^ 
12350 > 3761. Furthermore, we choose n = (p + l)/3. We have 

19 p _|_ x 2 
(a + «C YO ^ - < 3^ + ^ < p2 ' 

and hence, with the same notation as above, Equation (6.1) holds also in the current case. 
We have 3 ^ a 2 + b 2 ^ 4, a + b ^ < 2p, f = §, § < f ^ § < §, hence 



a2 + »2 a + & 



ci2 a 
T + 3p 



b_2 b 
3 3p 



= 1-1-0 = 0. 



3 3p 

Consequently, again, the p-adic valuation of j- ( an ^ >n ) equals —1 for our choice of p and 
n, which in particular means that 3^~[( an a ^' 1 ) is n °t an integer. 

Next let again (a, b) e {(2, 1), (2, 2)} + (3Z) 2 , but < b ^ |a. Since we have already 
verified the claim for a, b ^ 13000, we may assume a ^ 9200. (If a < 9200, then the above 
restriction imposes the bound b ^ |9200 < 13000.) Here, we choose a prime p = 2 strictly 
between |a and yjja. Such a prime is guaranteed to exist by Lemma 6.1, because, due 
to our assumption, we have |a ^ 7360 > 3761. Furthermore, we choose n = (p + l)/3. 
We have still 3 < a 2 + b 2 ^ 4, moreover 2p < f p < ±§a ^ a + 6 ^ fa < 3p, ^ = |, 

3 ^ 3« ^ 12 ^ 3' 3 ^ Ifln ^ Fin ^ 15p ^ 12 ^ 3' U-VLILV 



3p ^ 12 

a 2 + &2 



+ 



a + b 
~3p~ 



lOp ^ 3p 

a 2 



a 

3 + 3p 



12 

b_2 b 
3 3p 



1 + X (b 2 = 2) - 1 - X (&2 = 2) = 0, 



implying again that the p-adic valuation of j- ( an ^ n ) equals —1 for our choice of p and 
n, as desired. 



The next case we discuss is 



Since we 



(a, b) e {(2,1)} + (3Z) 2 and \a < b < 3a. 
have already verified the claim for a, b < 13000, we may assume b ^ 13000. Here, we 
choose a prime p = 2 strictly between -f^b and ^6. Such a prime is guaranteed to exist by 



Lemma 6.1, because, due to our assumption, we have ^b ^ 3900 > 3761. Furthermore, 



we choose n 
1 

Vr, 



3n 



(p + l)/3. Here we have 

an + 6n N 
an 





(a + 6)n 




an 




bn 


( 


p e 




pi ^ 




pt 



-1 + 



(a + b)n 




an 




bn 


pi 




pi 




pi 



+ 



a 2 



+ 



a + b 




a 2 


a 




b 2 b 
3 3p_ 


3p 




~3 ~ 


3p 
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In this case, due to the estimations 
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p 2 < 6| < (a + 6)n = (a + 6)^- < y&^y- < + 1) < 2p 2 , 



for p > 20, 



and 



we have 



T) 

p 2 <b- <bn < (a + < 2p 2 , for p > 20, 
3 



1-0-1 = 0. 

120 



(a + 6)n 




an 




bn 


p2 




p2 




p2 



(12 

3 



2 1 < A ^ Jl ^ 



Moreover, we have a 2 + b 2 = 3, 3p < b < a + b ^ —b < -^-p < 6p, ., ., , < . 

5b_ < 50 &2_1 1 / A < io th + 
21p ^ 63' 3 — 3> 1 ^ 3p ^ 9 ' bU LlldL 



a 2 + 62 O + & 

3 3p 



a 2 a 
T + 3p 



b_2 b 
3 3p 



2-1-1 = 0, 



implying again that the p-adic valuation of — (° n a * 6n ) equals —1 for our choice of p and 
n, as desired. 

The last case to be discussed is (a, b) G {(2, 1)} + (3Z) 2 and 3a < 6. Again, since we 
have already verified the claim for a, b ^ 13000, we may assume b ^ 13000. Here, we 
choose a prime p = 2 strictly between ^-6 and |6. Such a prime is guaranteed to exist 
by Lemma 6.1, because, due to our assumption, we have > 3761. Furthermore, we 
choose n — {p + l)/3. Since 

4 p + 1 44 
(a + 6)n < -&— ^— < ttP(p + 1) < P 2 , for p > 44, 
3 3 45 

for the p-adic valuation of 3~r( an a ^ n ) there holds again (6.1). In the current case, we 
have furthermore a 2 + b 2 = 3, a + b ^ |6 < ||p < 3p, ^ 



2/6 

3 



3' 3p ^ 9p ^ 45 ^ 3' 3 3' 



< tt ^ I!-, and hence 

3p ^ 15 ' 



a 2 + & 2 a + & 
3 + 3p 



a2 a 
T + 3p 



62 6 
3 3p 



1-0-1 = 0, 



implying also here that the p-adic valuation of g^— j- ( an ^ n ) equals —1 for our choice of p 
and n, as desired. 

This concludes the proof of the theorem. □ 



7. Concluding remarks and open problems 

In the proof of Theorem 2.1, assume that s is the smallest positive integer such that 
(a + b) | (p s — 1). Then we obtain the stronger inequality 



(7.1) 
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It is easily seen that s \ ip(a + b). However, such an upper bound is still likely much larger 
than the exact value of f(a,b) given by Sun [26]. For example, the inequality (7.1) gives 

7(7,36) <cl!^ = 2736, 

5 25 - 1 

/(10, 192) < = 1475362494440362, 

/(ll, 100) <;H^ = 15960, 
7(22,200)^^-^ = 7980, 

1000400 _ i 

/ (1999, 2011) ^ ss 5.272 x 10 1316 . 

J v ; 4010 

It seems that Theorem 2.2 can be further generalised in the following way. 

Conjecture 7.1. Let a and b be positive integers with a > b, and let a and (3 be integers. 
Furthermore, let p be a prime such that gcd(p, a) = 1. Then for each r — 0, 1, . . . ,p — 1, 
there are infinitely many positive integers n such that 

'an + a 



, = r (mod p). 
bn+ ° • 

In relation to Theorem 2.3, we propose the following two conjectures, the first one 
generalising Theorem 2.4. 

Conjecture 7.2. For any odd prime p, there are no positive integers a > b such that 

'an* 



. — (mod pn — 1) 
bnj 

for all n ^ 1 . 

Conjecture 7.3. For any positive integer m, there are positive integers a and b such that 
am > b and 

/ amn\ 
. =0 (mod an — 1) 
bn J 

for all n ^ 1 . 

Note that the congruences (2.1)-(2.6) imply that Conjecture 7.3 is true for 1 ^ m ^ 5. 

It seems clear that, for each specific prime p, a proof of Conjecture 7.2 in the style 
of the proof of Theorem 2.4 in Section 6 can be given. On the other hand, a proof for 
arbitrary p will likely require a new idea. 

We end the paper with the following conjecture, strengthening the last part of Theo- 
rem 3.1. 

Conjecture 7.4. For all positive integers n and non-negative integers k with ^ k ^ 
125n 2 — 25n + 4, the coefficient of z k of the polynomial 



[1 -gl0n-l)(l_ g 15n-l) 



30n 
5n 
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has non-negative coefficients, except for k = 1 and k = 125n 2 — 25n + 3 ; in which case 
the corresponding coefficient equals — 1. 
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